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We study the decay modes of a new, light spin-0 particle, arguing that if the mass of the 
(pseudo) scalar is ~ 11 — 15 GeV, it can have an appreciable branching ratio into bottomo- 
nium, in particular the rare rj^s. Using non-relativistic QCD (NRQCD), we calculate its 
decay rate to bottomonia for mass splittings greater than the typical momentum transfer 
within the bound state. It can exceed that of decays to other Standard Model fermions 
under the assumption of couplings proportional to those of the Standard Model. At smaller 
splittings, where our computational methods break down, we estimate the rate into bot- 
tomonia using data-driven methods. When the spin-0 state decays to bottomonia whose 
mass is too light to produce 5-meson pairs, we get a qualitatively new experimental signa- 
ture, decays to 6-quarks invisible to 6-tagging. Such a light, spinless particle can arise in 
extended Higgs sectors, making this channel potentially observable in decay chains initiated 
by the subdominant decay of a Standard Model-like Higgs to a pair of them. 



I. INTRODUCTION 



The LHC will soon cover the entire light Standard Model (SM) Higgs region. There is 
already convincing evidence that both Atlas and CMS are standing on the edge of discov- 
ery of a Higgs-like particle with mass close to 125 GeV jl, 2|. However, even this excess 
is confirmed as a new particle, many fundamental questions about electroweak symmetry 
breaking will remain unanswered. It will then be important to perform precision measure- 
ments, which would confirm or reject the hypothesis of the SM Higgs, probing various decay 
channels and the total width. Another crucial question to address will be whether it is a 
standalone particle, or a part of much bigger Higgs sector, as predicted by various beyond 
the SM (BSM) models, e.g. SUSY, or Higgs as a pseudo-Goldstone boson (pGB). 

Certain extended Higgs sectors allow for particles significantly below the weak-scale. Usu- 
ally, it is a pGB of a new, approximate global symmetry. One can find an illustrative example 
in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), which has an approx- 
imate R-symmetry in the limit of vanishing gaugino masses and A-terms. Therefore, if the 
NMSSM has small A-terms, it naturally contains a pGB of the continuous R-symmetry jsj. 
Additionally, little Higgs theories provide an example where the entire Higgs sector arises 
as the pseudo-Goldstones of a spontaneously broken global symmetry (see |4| for an early 
incarnation and |5| for a review and further references). 

We take this new, light particle, a, as part of the Higgs sector, and therefore expect 
its couplings to be roughly proportional to the SM Yukawas. Therefore, if m a > 2mb, 
its dominant decay mode would likely be a — > bb, while for m a < 2m>b, it will go mostly 
via a — > t + t~ (however, more exotic scenarios have been considered, e.g. a — >• cc 6| and 
a — > gg [7]). In this paper, we emphasize another non-standard decay mode, which was 
largely neglected in previous studies (see however a discussion in [S]). We consider m a large 
enough to evade the constraints from decays of bottomonia at CLEO and BaBar (m a > 9 



GeV [9!]), but not too heavy, such that direct decays of the (pseudo) scalar into bottomonia 
might still have appreciable branching fraction. As we will see, this corresponds to the 
mass region around ~ 11 — 15 GeV. We present a scenario where an important decay 
mode of a (pseudo) scalar is a — >■ bottomonium + X, where the S- wave pseudoscalars ?7h(n), 
3 Pj states Xbj( n ), or 1 D 2 (n) are the principle bottomonia involved (n denotes the various 
radial eigenstates) and X are light hadrons. We also show, that if a is a scalar, rather 
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than a pseudo scalar, this rate is suppressed. This suppression is less than one order of 

magnitude and therefore may be difficult to observe. Nonetheless, a discovery of the mode 

a — > bottomonium + X might serve as a hint of the parity of a. 1 

The collider physics of, say, a — > r]b + X decays is challenging, but interesting. First of all, 

the a can appear through non-standard decays of the SM-like Higgs boson. Theoretically 

it is even possible to achieve a spectrum, where h — >• aa is a dominant Higgs decay mode, 

getting the so-called "hidden Higgs" (see HQ for review and references therein). Given 

very strong hints for a SM-like Higgs at 125 GeV, this scenario is somewhat disfavored. 

Nonetheless, it is still possible that the process h — >■ aa is present, but subdominant to the 

standard h — >■ bb and h — > WW* decays. In fact, if the observed excess at 125 GeV turns 

out be a SM-like Higgs, then the gluon fusion cross section in the SM is ~20 pb, and so with 

20 fb _1 of data, one would have 400,000 Higgses. An 0(1 — 10%) branching ratio of h aa 

would therefore contain thousands of BSM events. If one were to observe additional heavy 

particles from the Higgs sector, other cascade channels would also open. For example, in 

I 

the NMSSM, one can produce the a in abundance thr oug h heavy charged Higgs decays [14j 



or in association with charginos and neutralinos [151. 116 . Alternatively, one can also look 
for as radiated from heavy flavors in bb and ti events 17j]. 2 

Quarkonia decay promptly and usually without leptons in the final state. Therefore, for 
those bottomonia that are below the two-S-meson threshold, the a ends up in two jets which 
cannot be 6-tagged. If the a comes from a rare Higgs decay, these jets can be collimated, 
resembling the situation in the "buried Higgs" [3]. Although the bottomonium decay mode 
is challenging, it has handles that make it different from a regular QCD jet, which we will 
elaborate in Section HVB1 The detailed analysis of rare Higgs decays in this channel, though, 
is beyond the scope of this project. 

The layout of the paper is as follows. In the next section, we give an introduction 



1 For current constraints on a light spin zero-particle see [s - 11 1 . We see that this possibility is wide open 
in the relevant mass range if the Standard Model-like Higgs is not too light. 

2 An alternative proposal to look for signs of the a is to study the r]b masses in detail to look for deviations 
due to mixing between the two. Additionally, one can look at the r-branching ratios of the T, attributing 
discrepancies with flavor universality to radiative T — > rjb'y transitions and subsequent 775 decays 18 - 2o| . 



We will have much to say below about the interaction between the a and the 775. However, given the 
observational rareness of the 775 and the theoretical uncertainties on their properties, it may well be that 
the a will allow us to measure properties of the % rather than the other way around. 
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to NRQCD and how we can use it to calculate bottomonium production rates. We also 
discuss a range of validity of this technique and show that the decay a — > r] b + X is non- 
negligible in this regime. In Section IHIl we discuss small mass splittings between the a 
and the bottomonium, where our NRQCD calculation is invalid. We estimate the order of 
magnitude of the a's decay rate to bound states through a particular, computable channel 
involving a — i]t, mixing. We also incorporate mixing into a computation of the rate for 
a — > gg over our range of interest. Readers wishing to skip to the punchline can look at 
Figs. (j3J) and (JSJ) for a final plot of decay modes and their branching ratios, showing where 
a — > bottomonium+X can be important. Finally, in the fourth section we discuss our results 
and conclude. A discussion of open-flavor decays to bs in the mixing regime is relegated to 
the Appendix. 



II. NRQCD CALCULATION 

For mass splittings between the light spin-0 particle, a, and the bottomonium greater than 
the typical bound quark momentum, (p& ~ rribV, where v is the quark velocity in the hadronic 
rest frame) we can compute the decay a — > bottomonium + X in NRQCD simply. Due to the 
factorization of long from short-distance physics, the calculation combines a perturbative 
QCD process with a nonperturbative factor taken from data or estimated through its power 
counting. After giving an overview of the effective field theory, we present each part of the 
computation in turn. 



A. NRQCD basics and validity range 

Non-relativistic QCD (NRQCD) is a rich example of an effective field theory where certain 
particles (c and 6-quarks) are light enough to exist as propagating degrees of freedom, but 
heavy enough that we can integrate out their pair creation, decoupling quark from antiquark, 
and can perturbatively expand in powers of their hadronic rest-frame velocity, v, as well as 



their inverse mass, 1/M, (we follow the treatment in Ref. 21]). Thus, we can work in the 



non-relativistic limit for the heavy quark fields, giving the following lagrangian, 

£nRQCD = £heavy + flight + SC, (1) 
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where £u g ht is the usual lagrangian for gluons and light quarks, 5C contains the correction 
terms that systematically give back full QCD, and 

= if (iD, + ^)i> + x* (iA - j£) X- (2) 

We have separated particle from antiparticle, and our fermions have been reduced to two- 
component Pauli spinors, with ip annihilating heavy quarks and x creating heavy antiquarks. 
The main advantages to using NRQCD for computing decays to bottomonia are two-fold: 
factorization and power counting. In fact, as we will show, we will need very little of the 
full NRQCD machinery once we understand how to make use of these two properties. 

The entire effective field theory approach is applicable because we have three parametri- 
cally separated scales in the problem: 

• The heavy quark scale, M. 

• The scale of momentum transfer inside the quarkonium, Mv . The size of the quarko- 
nium is also characterized by this scale, r ~ (Mv) . 

• The quark kinetic energy scale Mv 2 , which is also an energy splitting between the 
radial excitations of quarkonium. 



For 6-quarks the typical velocity of the quark in the bound state is v 0.3 [22 1. 

The parametric separation between the quark mass and the bottomonium energy and 
momentum scales lets us factorize its production into a short-distance computation, where 
we can work with perturbative QCD, and a long-distance, nonperturbative part, which takes 
the form of an expectation value of an NRQCD operator and accounts for the binding into 
a hadron. 3 For the case of a particle, a, decaying into bottomonium with a modest mass 
splitting between them, the leading contribution is: 

T[a^H + X] =J2na^bb(n,8)+g}(0%), (3) 



The validity of factorization for the production of quarkonium in NRQCD has never been proven and there 



is some evidence that the nonperturbative matrix elements are not strictly universal 23|,l2J|- Nonetheless, 
NRQCD has had quantitative success in calculating production rates (e.g. [3]). Since we will be taking 
estimates for our matrix elements of interest anyway, these concerns are beyond the order to which we 
are working. 
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where n labels the angular momentum ( 2S+1 Lj) of the bb pair, 8 denotes the color state of the 
bb pair and X is whatever hadrons emerge from the gluon emission as well as the conversion 
of the color octet bb state into a singlet bottomonium. Additional contributions, which 
are higher order in perturbative a s and have 3-body suppression, but can get kinematic 
enhancements can be important for large mass splittings (see 26] for such an effect in Z 
decays). The T term is a straightforward partonic calculation, but we need some way to 
determine (O^), which includes all nonperturbative effects. In some cases, (O^) are directly 
extracted from experiment, but not all matrix elements have been measured. In fact, none 
of those relevant for the a decay are known. Therefore, in Section IHBl we will use simple 
power-counting and comparison to measured charmonium rates to estimate those of interest. 

Before proceeding to the details of our NRQCD calculation, let us comment on its range 
of validity. Our perturbative calculation of F in Eq. [3] assumes that one emits a gluon 
harder than those exchanged between the bound quarks. Otherwise, one no longer has the 
separation in scales necessary to factorize the effects of binding into an overall multiplicative 
factor (cf. Fig. [TJ). Momentum of the softer, "potential" gluons is 0(mbv). Therefore, we 





FIG. 1: The decay of an a to 6-quarks with mass splitting < 0{mbv). In addition to emitting 
radiation, the bs exchange gluons that will bind them into a bottomonium. In this region of phase 
space, we cannot completely separate off the binding into the nonperturbative matrix elements of a 
small number of operators, but must consider the exchange of multiple gluons coming from infinite 
towers of operators. 

should demand 

n^a ^onium 

~ p(g) > m b v . (4) 
Numerically, m^v ~ 1.4 GeV and so we use the NRQCD rate given by Eq. [3] for splittings 
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greater than around this value. Below this scale, the parametric expansion in powers of v 
breaks down, and one needs to sum infinite towers of operators in order to compute. 

B. Calculation 

We take as our starting point an augmented Higgs sector containing a light (pseudo) scalar, 
a, with m a > 11 GeV. This means its splitting with respect to the lightest bottomonium, the 
77^(1) , is greater than 1.4 GeV, where we begin to trust the results of perturbative NRQCD 
(c/. Eq. H]). We now consider the perturbative and nonperturbative contributions to the a 
decay in NRQCD in turn. 

1. Perturbative portion 

We can have scalar and pseudoscalar couplings of a to fe-quarks, 

C-alb = Vbabb, or (5) 
C-alb = iybab^b , 

and similarly for other SM fermions. The origin of a in the Higgs sector motivates its 
having a large coupling to the third generation. Strictly speaking, we assume that the a's 
interactions are controlled by the SM Yukawas, and therefore other modes, like cc, t + t~ 
are also allowed. We take these channels into account in our rate and branching ratio plots 
(Figs. H]and|5]). 4 The large coupling of third generation fermions gives us, along with the 
usual, open-flavor bb decays, a — > bottomonium + X. This latter decay mode is always 
present, even if it is highly suppressed. However, it can become an important subleading 
process if the m a becomes close to the mass threshold of two B mesons. 

We first calculate the perturbative decay rates of a (pseudo) scalar a into the octet bb 
state and one gluon (see Eq. [3] and Fig. [2]). Despite the velocity suppression that arises 
from projecting the octet quark combination into the singlet, physical hadron, this is the 

4 Our use of these couplings is not an endorsement of a specific model, but rather an attempt to determine 
at the 0(1) level how branching ratios from the Higgs sector could appear. Specific scenarios {e.g. large 
tan(/3) NMSSM reducing cc) can suppress or enhance the rates we depict. However, generic scenarios 
minimally affect the relative ratios between them. 
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2S+1 Lj(8) 



+ 



<i 




2S+1 Lj(8) 



FIG. 2: Leading order diagrams for a —> bottomonia + X. Subsequent nonperturbative gluon 
emissions neutralize the color in both the perturbatively emitted gluon as well as the color octet bb 
pair. Unlike the situation depicted in Fig. Q3 the effects of binding factorize from those of radiation, 
and the former are captured numerically by expectation values of NRQCD operators, which appear 
as multiplicative factors (cf . Eq. [3]) . 

dominant contribution for decays of a light (psuedo)scalar. 5 After projecting the bb octets 
onto states with well denned angular momentum, we find that the (L = 0, S — 0) and 
(L = 1, S = 1) modes vanish. One can explain this behavior with simple C-invariance. 
This may be somewhat surprising as neither the color-octet bb pair nor the gluon is a C- 
eigenstate. However, the leading order diagram (Fig. [2]) is identical to a process in QED, 
with color factorized off as an overall factor, and thus we have an accidental C-symmetry at 
tree level. Therefore, only those processes that would be allowed for bb singlets and photons 
will occur. We determine the final state C-eigenvalue as if the constituents were themselves 
singlets, with each gluon bringing a factor of —1, and each bb pair a (— 1) L+S . Since any 
spin-0 particle is C-even, this leads to a selection rule 



thus ruling out 1 Sq and 3 Pj bb states as perturbative decay products. 

The nonvanishing modes have non-uniform dependence on the mass of the (pseudo)scalar. 
Defining for simplicity 



5 If in a mild abuse of notation we refer to this leading process as a — > bottomonium +g, then other 
important processes are a — > bottomonium +bb, a — > bottomonium + gg, and a — > bottomonium + qq, 
where q are light quarks, if we have sufficient phase space for the final state. For the decay of the Z into 
charmonium, the ipqq and ipcc are actually dominant, receiving a kinematic enhancement [3]. 



L + S = odd , 



(6) 



onium 



(7) 
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we get 

T(scalar) oc 1 — £ (8) 

r(pseudo, L = 0, S = 1) oc 1 - £ (9) 

r(pseudo,L= 1,^ = 0) oc \_I (10) 

One thing to note is that the 1 Pi pseudoscalar decay channel grows as m a approaches 
m onium while the other channels vanish. The origin of this structure is ultimately in the 
soft singularity of gauge boson emission. As we are expanding our bb pair in "spectroscopic" 
states, only those that do not require relative linear or angular momentum on the part of the 
gluon can get enhanced. The former constraint rules out having a pole for S'-wave bb states, 
since we require the quark pair to have identical momenta, and thus the gluon must recoil 
with finite momentum. The angular momentum constraint lets us see that enhancement for 
P-wave can only occur for a pseudoscalar a. The P-wave state allowed by C-invariance is 
1 Pi. For a three-body system, a given orbital angular momentum eigenstate has parity 

iVbody = PlP 2 P 3 (-l)'(-l) L , (11) 

where in our case i is the orbital angular momentum of the bb pair, and L that of the gluon 
with respect to the center of mass of the pair. The P, are the intrinsic parities of the decay 
products. A pseudoscalar allows L — 0, while a scalar does not. It may concern the reader 
that we have a pole in our formula associated with an IR divergence. However, in the regions 
of phase space where we use NRQCD, Am > 1.4 GeV, and the (1 — £) _1 factor only gives 
a modest enhancement, ~ 3 — 4. Long before the decay rate can blow up due to this term, 
one must account for both the exchange of "potential" gluons that bind the fermions as well 
as the reduction of possible hadronic final states from Aqcd- We estimate the decay rates 
for mass splittings which must take account of these effects in Section IIHI and find that the 
rate indeed turns over instead of blowing up. 

At last we present the explicit results of our calculations. The decay rate that we find 
for the pseudoscalar is 

I> -> i + X) = ^4 ( (1 - "hW*i)> + 4^^(^( 1 P 1 )) N ) , (12) 

where i stands for an arbitrary bottomonium state, which we will further specify in the next 
subsection. We will also discuss systematically why these matrix elements are the leading 
contributions in the m&, f-expansions. 
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The analogous result for a P-even scalar reads 

r(a P _ even -> i + X) = (1 - [m? (^( 3 ^i)) + 4 (OlCPi))} ■ (13) 

Even though this result is suppressed compared to the pseudoscalar, it is nonetheless in- 
teresting as it is part of the bottomonium decay rate of the Standard Model Higgs particle. 
Just as [26) found with the Z though, for a 125 GeV Higgs, it will be subdominant by a 
couple orders of magnitude to processes where a b or light quark fragments into a bottomo- 
nium, h — > (bottomonium + bb) or h — > (bottomonium + qq). While difficult to observe at 
the LHC, these processes potentially be measured in future linear colliders, allowing us to 
quantify the nonperturbative physics involved. 



2. Nonperturbative portion 

Even though Eqs. f|T2|) and ffT3l) are exact at leading order in QCD and are the lowest-order 
contributions in the and f-expansions of NRQCD, it will be impossible to translate them 
to quantative predictions if we know nothing about the long-distance physics incorporated 
in the expectation values, (O). 

For pseudoscalar decay (Eq. [12]), the perturbative portion of the 1 P\ term is relatively 
large, so we ask which physical bottomonium has the biggest overlap with the |66( 1 Pi), 8). 
One example is pseudoscalars, namely f]b{n). As with any other quarkonium state, 1%) is 
a superposition of Fock states, including bb pairs and possibly gluons. Adding a gluon or 
flipping a spin suppresses the probability to find that state by additional powers of v. Here 
we expand |?7&): 

\Vb) = \bbCs )i) + 0{v) \bbCPihg) + O (v^ 2 ) \bb{ 3 S 1 ) 8 g) + O (v 2 ) \bbCs ) 8 /igg) 

+ 0(v 2 ) \bb( 1 D 2 ) s/1 gg) + ... , (14) 



using the power counting of 21|. 6 Each electric gluon in the state brings with it a factor 
of v , while the magnetic gluon (found in kets where the quarks are spin-flipped) has an 
additional suppression, for an overall factor of v 3 ^ 2 . 



(i 



This expansion closely follows a similar expansion for Z decaying into T 26| , but uses the slightly different 



power counting assignment of 21 1 
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The r)b thus includes the |66( 1 Pi), 8) Fock state at 0(v), which is the lowest possible order 
the physical singlet can mix with the octet fermion pair, including one dynamical gluon to 
get an overall singlet. Therefore, a — > T]b + X is a dominant bottomonium decay mode and 
we will be interested in the value of (02 b ( 1 Pi)) to estimate its rate. Unfortunately, r/f, is not a 
well observed particle and we cannot take the relevant matrix element from data. However, 
we can make use of the power counting of NRQCD to estimate the size of nonperturbative 
contributions with factors of m& and v. As a cross-check, we can also take measured values 
from the charmonium system and rescale them appropriately for bottomonium. 

We begin by understanding how it is that this matrix element describes the process of 
turning the pair of 6-quarks into an 775. The operator Ofi^Pi) creates a bb pair, computes 
its overlap with the asymptotic hadronic state (r]b + X\ and takes the outer product with the 
complex conjugate state so that its vacuum matrix element gives the probability to produce 
the 7]j, (plus whatever else) from the heavy quark pair. We have four fermion fields, which 
places this operator in 8C, in terms of Eq. [TJ This makes sense as the creation of bb pairs is 



short distance physics from the perspective of NRQCD. From [21], we have 

OffPi) = x ] (-^) {«) ^ (-^) ( 15 ) 

where 

{al b a Vb )=J2\Vb + X}(v b + X\. (16) 
x 

We are now in a position to estimate the size of the matrix element, taking the state and op 
erator normalization of 21 



271 ] . The size of its contribution is m\v 7 , the breakdown for which 
is given in Table [H Numerically, ((9g b ( 1 Pi)) ps 1 GeV 5 (hereafter we use m;, = 4.7 GeV). 
We see furthermore, why we have truncated Eqs. ffl2l and ffl3|) at two operators. Including 
higher orbital angular momentum states brings further v suppression, as the corresponding 
four-fermion operators must have additional covariant derivatives. 

Alternatively, while this matrix element has not been measured for bottomonium, a 
related matrix element in charmonium, (Og^( 3 Pj)) has been determined 28, 29], albeit 



in linear combination with (Og '( 1 So)) ■ Assuming the two terms are comparable, we have 
KT 2 GeV 5 w (O^^Pq)) ps 1(0^(^1)). 7 The second approximate equality follows from 



7 Recently, groups have attempted to determine (Og^( 3 Po)) 30-32J and even (Oj{ 3 P )) [33| , treating 
QCD at NLO. This is a highly nontrivial analysis on both the theoretical and experimental sides. Un- 
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Factor 


Origin 


(mbv)~ 3 


Volume factor from 
operator spatial integral 


(rribv) e 


4 heavy quark fields 


o 
V 2 


Overlap of 66( i J Pi) 8 

with % Fock state (c/. Eq. I14p 


(m b v) 2 


D l in operator 


x 7 

miv 


Total 



TABLE I: Power counting for (Of { l Pi)} 



spin symmetry, which holds that properties of quarkonia that differ only by flipping spins 
should be equal up to 0(v 2 ) corrections from chromomagnetic operators once the ratio of the 
number of spin states is accounted for (hence the factor of 1/3). For charmonium, v 2 = 0.3, 
and we can therefore convert from charmonium to bottomonium by taking 

<i 7 

mCPx)) « ^{Of{ 1 P 1 )) » 0.3 GeV 5 . (17) 

m c V c 

Thus, for our computations below, we will split the difference of our estimates and take 
(Of^Px)) = 0.5 GeV 5 . Taking this number (and the same procedure for other nonpertur- 
bative estimates) and the factorized equation ([3]), we can proceed to numerically estimate 
the decay rate from Eq. (JT2]) . which we plot in Fig. HI 

The 7]b is not the only particle which mixes with the |66( 1 Pi),8) state at the leading 
orders in v. We show the full list of these states in Table [TTJ Although hf, and 1 F 3 will also 
have non-vanishing overlap with this states, we will further neglect them in our discussion 
since technically they overlap at higher order in v. Nonetheless, the 1 D2 state should have 
the same overlap with |66( 1 P 1 ), 8) as rjb and must be considered as well. It is straightforward 

fortunately, these results have large uncertainties and currently take negative central values. Since the 
a decay rate is proportional to the nonperturbative matrix element, not its mod squared, we cannot use 
these determinations. Were these negative values to persist, the thought experiment of coupling a light 
pscudoscalar a to b quarks provides evidence for the breakdown of factorization in NRQCD, recently 
discussed in [24j]. The nonperturbative contributions to the a decay to bottomonia would have to be 
positive. We therefore take our estimates of the long-distance matrix elements from power counting in 
bottomonium and from the order of magnitude provided by the charmonium data, which are in rough 
agreement. 
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to see why this has the same power counting. Like the r)b, it has no relative spin-flip relative 
to 1 Pi, just a change in orbital angular momentum by one unit. 



Hadron 


0{v) 


Vb 


V 


X D 2 


V 


XbJ 


v 3/2 


h 


v 2 


'F 3 


V 2 



TABLE II: Hadrons whose Fock state expansion contains a 66( 1 Pi)s component, and the order at 
which it appears (c/. Eq. (THl) ). 



Similarly, we can understand the Xfej' s power counting as it has the same orbital angular 
momentum as 1 Pi, but differs by spin-flip as it is 3 Pj, bringing in an extra factor of v 1//2 in 
the amplitude. 8 Using the results from Table flU we get 



(18) 



The 3 in Eq. [18] comes from the suppression of the Xbj rate due to an additional power of 
v. However, since there are three nearly mass-degenerate Xbj states, it is effectively leading 
order as well. We base our numerical estimates for the other states on the rji,, with (Og 6 ( 1 P 1 )) 
= 0.5 GeV 5 . For the unmeasured r]b, Xbj and all 1 D 2 masses needed for Fig. H] and [5J we 
use the values determined using the potential model in Ref. j^. 9 Eqs. fT2l and [TBI involve an 
additional nonperturbative term, (0 8 ( 3 Si)). Of our states of interest, the dominant matrix 
elements are (0$ bJ ( 3 Si)) ~ 0.3 GeV 3 , with J = 0, 1,2, which we have estimated by power 
counting. 



The power counting of 2l| used in Table ITT1 assumes that itiqv 3> Aqcd- A more conservative approach 



takes ttiqv ~ Aqcd, as is done in 



34 



361] . These parametrics seem disfavored for bottomonium, which is 
why do not use them in our analysis, but may relevant for analyzing charmonium, since m c v w 700 MeV. 
For completeness though, we mention the effect adopting them would have on our results. In this regime, 
one no longer has a hierarchy between electric and magnetic gluons, and so subleading Fock states with 
a single magnetic gluon are only suppressed by 0(v) 36J. This means that each Xbj angular momentum 
state would be the same order as r\h or 1 Z?2 and collectively they would dominate. Furthermore, in the r\b 
itself, the 3 iSi component of its Fock state (c/. Eq. [T4"]) would come in at 0(v), the same as 1 P\. 



Practically we take from data only 775(1) resonance mass, however it has been recently claimed in 381 ] . 
that the 775(2) has been also experimentally observed. 
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Before proceeding, we make one final comment on the decay of scalar vs. pseudoscalar a. 
Since the pseudoscalar rate has a pole in the 1 Pi term for m onia /m a = 1, while the scalar 
rate vanishes in this limit, the former will be larger. For m a = 12 GeV, the ratio is a factor 
of 4. Should we find a spin-0 state in this mass range, improving on our tree-level, estimated 
rates, especially by an independent measurement of nonperturbative matrix elements, may 
allow a determination of the P eigenvalue of a. 



III. BINDING REGIME 



For mass splittings > 0(mf,v) between our light, spin-0 particle, a, and our bottomonium 
decay product, NRQCD allows us to easily compute the a's decay rates in Eqs. f[T2"j) and ([TBI . 
As a local effective field theory, we can use the expansions in a s , nib, and v to systematically 
capture leading contributions with a handful of operators. For smaller mass splittings, the 
calculation becomes more involved. We argued in Section III Al that we expect a breakdown 
when the momentum of the emitted gluon becomes comparable to the momentum transfer 
within the bound state. We can see this operationally as follows. In performing the NRQCD 
calculation, we expanded about the limit of zero relative quark momentum, which gave us a 
series in v 2 associated with higher orbital excitations. However, as we squeeze the sup por t 



of the amplitude into a region of width emj, our expansion switches to powers of v 2 /e 39]. 
Thus, for energy splittings of 0(mbV 2 ), which correspond to a mass splitting of 0(m&i>), we 
can only calculate if we can sum up infinite towers of operators built on those giving the 
leading contributions for larger splittings. Furthermore, we see that for even smaller mass 
splittings (e <C v 2 ), we lose calculational control entirely, as suppressed operators become 
important, since v n (v 2 /e) m > 1. 

We do not do so in this work, but the OPE provides a means to sum the infinite towers 
of operators necessary for Am ~ Oim^v) (for a review of OPEs in 6-physics, see j^j, and 



an OPE for bottomonium decay is found in 4l|]). 10 We define 

T = i J AT(4 5t (i)£ a5t (0)), (19) 



10 An alternative method for calculation could be to use the effective field theory potential-NRQCD (pN- 
RQCD). It is designed to operate at mass scales < 0(mbv), where one integrates out the gluons responsible 
for binding, those with (E, p) ~ (mt,v 2 , rribv), leading to a spatially non-local, instantaneous potential 



(see 42|,|43( for reviews). 



14 



where C a i b is one of the terms in Eq. |6j Then, by the optical theorem, 

(20) 



^S.states+x _ Im(a|T|a) 



m„ 



Performing the OPE on T gives us an infinite series of NRQCD operators whose expectation 



values we can sum into nonperturbative structure functions, as in [39l. |41|. 

Although we do not perform the OPE calculation (which is expected to give a rather 
precise answer for Am ~ 0(rribv)), we roughly estimate the decay rate to bottomonia using 
a different technique. We notice that when m a becomes comparable to or smaller than the 
heaviest rjb mass (=11 GeV), certain decay rates will receive significant effects from a — rj^ 
mixing. These are 

1. a -»■ gg 

2. a -> r] b + X 

3. a — > (open-flavor) bb 

We will spend the rest of the section discussing 1) and 2), while leaving 3) for an Appendix 



(since we use the approach of [44J, albeit with updated data). After we compute the mixed 
content of physical states, the a — )■ gg rate will follow readily. We compute rate a — > + X 
solely due to hadronic transitions between bottomonia after mixing with the a. This will 
give us an order of magnitude estimate on the rate for Am < 0{rribv), where the calculation 
is either involved or very difficult. Since we are in a region where parametric control is 
breaking down, it is difficult to quantify how much of the total rate these mixing effects 
represent. However, at those locations where the a is degenerate with an rjt, and mixing is 
maximal, we should expect the decay rate due to mixing to dominate, and thus these points 
give us an approximate upper bound. 



A. a — rib Mixing and gg annihilation 

Since the a has the same quantum numbers as the ?7feS, it will mix them. 11 We can calculate 
the mixing from the diagram in Fig. |3j Quantifying the conversion of a 1 Sq bb singlet into an 



Since parity is broken in the electroweak sector, the a will mix the the parity- even Xbo as well, but this 
rate is suppressed by weak coupling, a loop factor, and the W mass so we neglect it. 
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FIG. 3: Feynman diagram for a — r/b mixing. The blob represents the nonperturbative contribution 
of the radial wavefunction at the origin, i? % ( n )(0). 



?7fe hadron does not require the full machinery of NRQCD, but can be described simply via 



the hadronic radial wavefunction at the origin 



44j . 12 The mass matrix that controls a — rjb 



mixing thus contains the following off-diagonal entries: 



5m a-%(n) 



Vb 



An 



-m r 



Vb(n) \Rv b (n){V)\, 



(21) 



One can obtain the nonperturbative quantity |i2^ b (0)| from experiment. While there is little 
data for the t^s directly, by spin symmetry we expect the wavefunction values of the T to be 
the same up to suppressed chromomagnetic effects. We can therefore get a numerical value 
for Eq. Ell from the PDG 



46| via the following formula: 



|i?T(n)(0)| 2 = r(T(n)->e + e- 

Eqs. ( |2T|) and ( |22|) give us the mixing matrix 
/ 



9 m\ ^ 
4a 2 



M 



\ 



m %(l) ~ im Vb{l) T r, b (l) 








16 a s {m r ] 
3tt 








-1 



(22) 



m v b (6)- V!l >" 



l Vb(6) r %(6) J 

We include six 775 states, corresponding to the six observed T particles. Other than the 



(23) 



Vb 
45 



2) 



47 



37 



jtheir masses are unmeasured and have to be taken from calculation (Refs. 
49J). We give their values, along with |i? r?ft ( n )(0)| 2 in Table HTT1 The various nonzero 



widths, T, give us a non-Hermitian matrix, but it is still symmetric and diagonalizable 
via a complex rotation among the a and rj^n). For the premixed widths, we take T a to 



12 In this section, we switch to working with a Schrodinger potential model (see e.g. [45I]). The mixing, anni- 
hilation, and hadronic transition rates we need can be written very easily in terms of the bb wavefunction, 
and the numerical value extracted from experiment. One can recast the wavefunction in terms of NRQCD 



matrix elements, e.g. for the 7/f, radial wavefunction, R nb {r) — J jj-(0\xH~~ r /2)V'( r /2)|?7b) 21\ 
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m Vb ( n ) (GeV) 


\R Vb[n) (0)\ (GeV) 3 / 2 


1 


9.4 


2.7 


2 


10.0 


1.9 


3 


10.3 


1.6 


4 


10.6 


4.7 


5 


10.85 


4.7 


6 


11.0 


3.0 



TABLE III: Masses and wavefunctions at the origin for the ijj, states. 



include decays to cc, t + t~, gg, and in the regime of validity for NRQCD (> 11.5 GeV), 
the bottomonium + X decays determined in Section [Til For the decay rates of the rjb{n) 
themselves, we consider three classes of decays: annihilation to gg, hadronic transitions 
to lower lying r/f,, and for the r]f,(5) and r/b(6), decays to B mesons. 13 Once again, our 
approximation is to neglect any decay of the a to bottomonium that does not involve mixing 
for m a < 11 GeV. 

Computing the decay to gluons is straightforward. The ?7 fe 's rate is 



r(r} b (n) ->• gg) 



a s (m 



3 m 2 , \ 



\R 



Vb(n) 



(24) 



We also have the a's "own decay" a — > gg. For the one-loop contributions to this rate, we 
use the formulas of Refs. (3, 3' including ti, bb, and cc loops. Diagonalizing the matrix in 
Eq. (|23|) and rotating to mass eigenstate basis, we get for our physical a particle 

a P hy S . = c a a +) j c i T] b (i), (25) 



where the q are generically complex, and we can get nontrivial interference effects between 
different decay channels, since these add at the amplitude level, 



M(a phys . gg) = c a M(a gg) + y^ j c i M{r], 



(26) 



In fact, the changes to masses of the physical eigenstates from the naive unmixed values 
are negligible, but we see the change to T(a -> gg) is pronounced (c/. Fig. [4J. However, 



13 Unlike the T system, where the AS famously decays to two B mesons, the rjhi^S) is disallowed by CP and 
P symmetry from decaying to BB, and kinematically forbidden to go to higher mass, higher spin mesons. 
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compared to the other decay rates of the a, the branching ratio to gluons is still small. For 
particular models though, if the effective coupling to r + r~ gets extra suppressions beyond 
the Yukawa coupling, this rate could be important and is used in the buried Higgs Q. 

B. Bottomonium hadronic transitions 

In this subsection we estimate decays of the a to bottomonia via its mixing with the rjb 
states (Eq. Eli) , which undergo a subsequent hadronic transition. In the case of decays into 
gluons we could easily calculate the a's "native" decay for m a e [9.4, 11] GeV. When the 
decays proceed into bottomonia this computation involves an infinite number of operators, 
as discussed at the beginning of this Section. We will therefore estimate the a's decay rates 
to bottomonia in this region from the mixing contribution alone. If we were in possession 
of exhaustive data cataloging the decays of all six r^s, then this would involve little beyond 
diagonalizing the mass mixing matrix (Eq. 1231) . Since only two 77&S have been observed and 
their decay channels barely determined, we will need an alternative means to estimate their 
decay rates, which we will plug into the formula, 

T(a — » bottomonium + X) = |q| 2 T(r]b(i) — > bottomonium + X)9(m a — m Vb ^). (27) 

i 

In addition to the approximation of dropping a decays that do not involve mixing, we also 
sum the rates with identical final states incoherently as the a only mixes strongly with one 
bottomonium at a time, and we ignore the kinematic difference between m a and Tn^.u\ since 
in the mass range where mixing is strong, this will be a small effect. The ^-function though, 
enforces that kinematically forbidden processes cannot occur. 

Although we do not have data on the exclusive channels of their spin-flipped partners, 
the Ts, are much better observed. Thus, we will use their observed hadronic transition rates. 
For the low-lying (n = 1 — 3) states, we can justify using Ts in place of ?7f,s by the QCD 



multipole expansion (for details see 50|, [51| and references). This relates transition rates 



to the wavefunctions calculated in a Schrodinger potential model, and has had quantitative 



success in predicting decay rates 



52 



2|-l57l|. We will use the same approximation as in Section 
IIII Al that the rjb and T wavefunctions will agree up to 0(v 2 ), spin-suppressed effects. Thus, 
we can take the T transition data, which we give in Table HVl 14 



For the low-lying states where the multipole expansion is more trustworthy, we can calculate the self- 
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For the higher T excitations, the multipole expansion fails badly, by as much as two orders 
of magnitude 



58 



59]. Interestingly, the hadronic transition rates for the higher states, are 



relatively large, 0(1) MeV, agreeing with what we find in the nearby perturbative NRQCD 
region (c/. Fig. Hj). It has been suggested that one must take into account the effect on- 
shell 5-meson states can have on transition rates. Attempting to parametrize these, the 



authors of Ref. 



58 



59] could predict rates at the factor-of-two level. More recently, data 



on the hadronic transitions of T(5) have shown the presence of nearby narrow resonances 



responsible for the enhancement of its decay rate |60M62| . The transition proceeds in two 
steps, T(5) —> Zbn followed by — > T(n)n (n = 1 — 3), where Zb is the new, charged, spin-1 
particle, likely a molecular bound state of B mesons. By heavy quark spin-symmetry, an 
analogous molecular state with spin-0 is predicted to exist, Wbo, with comparable couplings 



to the r]b as Zb has to the T 63M65| . Thus, we should expect the same enhanced decay 



rates for the pseudoscalar bottomonia as for the vectors. This justifies using the T hadronic 
transition data for the higher excitations (which we present in Table HV|) in Eq. [27] for the rjb 
decays. 15 In detail, due to spin and kinematic factors there may be 0(1) differences in the 
rjb and T rates, but for the order of magnitude estimate we are trying to get for a decays in 
the binding regime, these are beyond the order to which we are working. 

For the higher states, there can also be sufficient phase space to decay through heavier 
pseudo-Goldstones such as 77 and K. We list these additional processes in Table HV] which 
can have comparable rates to the dipion modes, but we do not include them in our branching 
ratio calculations since we cannot do so systematically, lacking data for the T(5, 6) data 
for the single- 77 and/or KK channels, and we are only attempting an order of magnitude 
estimate. 

In the end, we find that the rate of a — > T]b(n)inT from hadronic transition after mixing is 
significant only where it has large overlap with the 775(5) and 775(6) states (cf. Figs. H]and|5]). 

coupling of the 77b (i)s to pions, ~ ??b(i)?7b(i)7r7r. Using the parametrization of 50] to go between gluons 
and pions, and the potential model of |54j |. we find, for example, the rate of a r)b(l)mr through the 
77b (1) self-coupling to pions ~ 10~ 8 GeV, an order of magnitude below any other rate we calculate, and 
far below the leading processes. We thus neglect these self-coupling mediated decays. 
15 Our use of data breaks down for the T(6), whose dipion transitions to other T states are unmeasured. 
Thus, for its rates we use a combination of those provided in 58] and an assumption that BR[T(5) — ► 
T(5 - n)\ « BR[T(6) -> T(6 - n)} for the others. 
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Initial state i= 


Ej r(T(i) -> T(j)^r) MeV 


Multipole valid 


Other important transitions 


2 


9 xl0~ 3 


Yes 




3 


2 xl(T 3 


Yes 


T(2) 77 


4 


5 xl0~ 3 


No 


1(1)7? 


5 


1.5 


No 


T(l)iOf, T(l, 2)t?? 


6 


3 


No 


T(l, 2)M?, T(l -3)r/? 



TABLE IV: Total hadronic transition rates for the dipion decays, T(i) — > T(j)7T7r, which we use 
for the 77;, transition rates in Eq. [271 

IV. DISCUSSION 
A. Rate Plots 



We present our results for the a's decay rates and branching ratios for the region 9-15 
GeV. For these figures, we restrict to the case of pseudoscalar a. Most simple to describe 
are those to open fermion pairs, partonic-66, cc, and t + t~. At tree level, this rate is given 
by 

' 1 (28) 



G) 



r(a->//) = N c ^ r m)m a [ 1- 



2v/2 



7T 



For the partonic calculation of quark rates, we also apply the NLO QCD correction given in 
Ref. 44j. For a — > gg, we use the result of Section fill Al where took into account the effects 
of a — rjj, mixing. Mixing also affects the turnover to production of 6-flavored mesons, where 
we have used a phenomenological interpolation (cf. Eq. IA4I) to go between a partonic QCD 
calculation for 6-quarks with m& = mg (Eq. |2"8|) and the hadronic mixing formula, which 
uses Eqs. IA1I and IA4I As discussed at the beginning of Section DTD and in Section IIII Bl for 
m a < 11 GeV, we estimate its decay rate to bottomonia only through its mixing with rjb, 
which undergo hadronic transition decays to give a — > ri^-n-n . We use the T transition rates 
as approximately those of the 77^ by spin symmetry in Eq. [271 Lastly, we have the process 
a — > bottomonium +X calculated systematically in Sec. Ill Bl for m a > 11 GeV, using Eq. [3j 
We see that the hadronic transition after mixing and NRQCD computations come within 
an order of magnitude where we switch from one method to the other. This gives us reason 
to believe a detailed study could find a smooth interpolation of the total a — > onia + X 
rate as one passes from the regime of Am a _ onia > m b v to Am a _ onia ~ m b v at values roughly 
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comparable to those we depict. 

In Fig. HJ for visual simplicity, we only show the rjf, NRQCD-derived decay rates, even 
though l D% and Xbj are large as well. For each of the six states, we cut off the curve when the 
a — rjb mass splitting becomes smaller than 0(rribv). To plot branching ratios in Fig. EJ we 
need a full accounting of the various decay rates, and thus include all of our leading NRQCD- 
derived rates. Since the 77^,(5, 6) decay overwhelmingly to 5-mesons, which experimentally 
look like bb decays, we do not add them, nor any of the other above-threshold bottomonia, 
to the "onia + X" curve. The NLO QCD computation we use for the decay to bb is really an 
inclusive rate for all hadronic channels that contain a bb pair and a perturbative gluon. This 
therefore includes the a — > bottomonium + X final state, and so for Fig. |5j we subtract our 
rate obtained from Eq. [TJ]for the sub-threshold bottomonia from the inclusive one (Eq. IA4|) 
to get branching ratios for open-flavor channels versus those that cannot give two S-mesons. 

For the sub-threshold bottomonia, we use the NRQCD calculation down to a mass of 
11.6 GeV, as this leaves a large enough splitting to the heaviest such state, r/b(4), that X 
can involve many possible decay channels. We have left the range [11.1, 11.6] GeV blank in 
Fig. El as we lack a compelling way to estimate our way through it. Since the mass splitting 
in this range is 0(rribv), it is possible in principle to use an OPE based in Eq. [19] to compute 
the rate, provided one can numerically determine the resulting nonperturbative structure 
function. It would be an interesting exercise to see how extending the NRQCD results down 
to this mass splitting compares with the black curve we have estimated from a — rjb mixing 
followed by hadronic transition. 

B. Observability 

Looking at Fig. |5]we see that the bottomonium channel is almost always subdominant, 
either to bb or to t + t~ and cc. However, in some parts of parameter space it is not negligible 
and can account for almost 30% of the entire branching ratio. Bottomonium does not look 
very different from a regular QCD jet (even though it typically has fewer tracks) and therefore 
one can view the enhancement of the bottomonium decay mode as an effective enhancement 
of the gg channel. Indeed, in the a rest frame, r]b is recoiling against a gluon jet, while rjb is 
probably impossible to distinguish from a gluon jet on event-by-event basis. Practically, b- 
tagging does not work on the hadronic decays of sub-threshold quarkonia as they annihilate 
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Decay Rate (GeV) 



0.001 




10 
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FIG. 4: The total decay rates for pseudoscalar a for each channel in GeV as a function of m a . The 
bb rate is taken from the interpolation equation (jA4j). For visual simplicity, we have only included 
the rib. We plot the hadronic transition rates from a — % mixing up to a mass of 11.1 GeV, at 
which point we switch to our NRQCD calculation, Eq. [T2l for the decay a —> %(1) + X. We do 
not use NRQCD for m a — m^ b t n \ < 1.7 GeV for any n. The rjb^Tnr curve is dashed as we have 
only included those contributions from a — rjb mixing. It is thus meant as an order of magnitude 
approximation, rather than a systematic calculation. 

rapidly to gluons (cf. [24]), and this is an emergence of "buried Higgs"-type signatures Q, 
but at masses far above the ditau threshold. 

Without further evidence for new physics in this particular window, it is perhaps 
premature to launch a dedicated collider study into the possibility of observing a — > 
bottomonium + X decays. However, there is still a benefit at this stage in a qualitative 
discussion of how this scenario would appear at the LHC, and how we might go about 
observing it. If we take the Higgs mass, rrih, to be around 125 GeV and m a = 12 GeV, 
then even in the Higgs rest frame the decay products of the pseudoscalar will be collimated, 
having AR ~ 0.3. For the region where a — > bottomonium + X decays are important, the 
gluon is fairly soft, having p < 5 GeV in the a rest-frame. 

A standard opening move for analyzing a jetty Higgs is to look at channels where the 
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Branching Ratio 
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FIG. 5: Branching ratios to the various decay products of an a pseudoscalar. We have summed over 
the four rjb , nine Xb J, and two 1 D2 sub-threshold states to get an overall a — > onia +X. We cease to 
consider a bottomonium decay of the a in NRQCD below m a = 11.6 GeV. We show the dominant 
contribution for smaller Am though, from the decays a —> rjbinr, starting at m a = 11.1 GeV, where 
we begin to see mixing effects with the heaviest r/&. The rjbjinr curve is dashed as we have only 
included those contributions from mixing. It is thus meant as an order of magnitude approximation. 
For the blue curve, we have subtracted the contributions of sub-threshold bottomonia from the 
inclusive bb + X rate given by Eq.Elto give the rate to open-flavored mesons. 



Higgs itself is boosted, coming from production processes such as V + h and tth 



then attempts to build one large "Higgs jet" either through fat cones [67j or the Cambridge 



66j. One 



Aachen algorithm 
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Having t 



his 



pileup and underlying event [66 



69 



et, and after performing some cleaning on it to remove 



701 ]. one can then look for structure consistent with 



having two dominant clusters of energy with approximately equal masses. For our case, just 



as in Ref. 



67j, the particle associated with the subjet, the a, is less massive relative to its pt 



than a QCD jet. Additionally, as in Ref. [67], we can exploit the absence of colored states 
until the a decays, which gives a different radiation pattern from prompt diiet production. 

Despite the similarities in technique appropriate to our scenario and [67]], and the simi- 
larity in final states, we have a couple additional handles. Firstly, in addition to having the 
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subjets at comparable light masses, their masses are near those of known states, rjb, XbJj 



and 
(of. 



37 



45 



ew of these have been measured directly, but calculations exist for the others 
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491 ]) likely to be valid to within 1 GeV, which puts them well within exper- 



imental resolution. Additionally, we can make use of the rj^s decay into a small number 
of hadrons. While the dominant decay products of the rjb have not been measured, the rj c , 
which similarly annihilates into gluon pairs, typically goes to just 2-4 particles. Thus, count- 
ing the number of tracks in the decay, and looking for a small number relative to a QCD 
jet with similar could provide useful additional information. Should it become necessary, 
one could even attempt to study in more detail the decay properties of these bottomonia 
states by combining NRQCD with the lattice. 

Lastly, we note that for many different arrangements of the a's coupling to SM fermions, 
over much of the m a range where bottomonium decays can be important, there are other, 
experimentally distinctive decay channels, one can observe. Below the open b threshold, one 
could search for events with h — >■ (r] b + X){t + t~). Above it, one could use 6-tagging on an 
a that decays to B mesons, and look for the other to decay to a light, sparsely populated 
subjet. Thus, despite putting the Higgs in the realm of jet physics, one has many handles 
for digging out the structure inherent to a decay chain passing through quarkonia. As a 
final aside, we mention a side benefit of having a scenario such as this realized in nature. 
Despite many decades of experimental S-physics, the rj b states remain very poorly studied 
despite being kinematically the simplest (spinless and S'-wave). Nature would have a wry 
sense of humor to deliver us these particles, and the other previously-unobserved mesons, 
out of its mechanism for electroweak symmetry breaking. 
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Appendix A: Open Flavor Decays 



The decay to open ^-flavored hadrons will dominate for sufficiently large m a . One must 
still determine though, where exactly in the m a range this process takes over. This is a 
somewhat involved question as it brings in issues of quark-hadron duality, IR divergences, 
and large threshold logarithms; we do not attempt to answer it in full. It is nonetheless 
interesting from the point of view of decays to bottomonia, as these can be an important 
subleading effect above the open-flavor threshold, and are kinematically allowed even where 
B meson decays are forbidden. The zeroth order approximation for a — > bb is to simply use 



the partonic formula with = m B , where B is the lightest B meson. Following Ref. |44|, 
we wish to improve on this naive treatment by incorporating effects due to a — 775 mixing, 
discussed in Section IIII Al in the threshold region for B meson production. We know the 
partonic calculation is making a kinematic mistake in the threshold region, and this is our 
attempt to partially correct it. The pseudoscalar decay into b quarks is an S'-wave process 
by parity. However, the decay into the lightest allowed physical mesons, a — > B*B + BB* 
is P-wave. Thus, we include the modification of the threshold rate due to the presence of 
f]b{5) and ?7&(6), with their appropriate P-wave coupling to P-mesons. 

One important qualitative difference between S and P-wave is the momentum dependence 
near threshold. The physical, P-wave, decay will be proportional to the decay products' 
momenta, which goes to zero at threshold, while the partonic S'-wave channel receives no 
such suppression. We write out the decays of rjb and T to the lightest B mesons so that we 
might use data from the latter to estimate the former, 

M( Vb ^B*B) = a B k5 X0 d° 0X (Al) 

M( Vb ^B*B*) = a B k5 x - x d° QXf 

M(T -> BB) = a B kd\. 

M(T -> B*B) = a B k\ f d\. Xf 

M(T ->■ B*B*) = a B kd\ iXf , 

and similarly for B — > B s . The helicities are A, for the decaying meson, A, A for the final 
state particles, and A/ = A — A. The norm of the three-momentum of the decay products 
in the rji, or T rest frame is k. This factor arises from the P-wave nature of the decays. 
To avoid introducing extra notation, we have used the same symbol throughout, but k will 
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differ with each process. The terms are Wigner's little-d functions. The overall factor 
as is a three- meson coupling and is assumed to be the same for all of Eq. I All since each line 
can be obtained from another via spin flips. 

In practice, the assumption that these different channels are all controlled by a single 
coupling holds only to within factors of two in the B decays of the T, and at the order of 
magnitude level for B s . We nonetheless make it as we can only hope to understand threshold 
effects up to 0(1) factors anyway. Additionally, the branching ratio of T(5S f ) to B mesons is 
nearly three times that to B s , so we accept a larger uncertainty in the subdominant process. 
We fix the couplings as and clb s such that the overall rate of T(5) into the three channels of 
Eq. IA1I agrees with data. The inclusion of B s decays is a change we have made from [441 ]. as 
experimental evidence for such processes in T(5S) decays was just starting to be reported 
at the time. Calculating the ratio of rates from Eq. I All we get 

r(?fo(5) -»■ B mesons) « 0.9T(T(5) -»■ B mesons) (A2) 

r(r/b(5) — > B s mesons) ~ 0.65T(T(5) — > B s mesons) 

r(77fe(6) —> B mesons) ~ r(T(6) — > B mesons) 

r(r/ft(6) — > B s mesons) ~ T(T(6) — > B s mesons), 

where we have summed the first two lines of Eq. IA1I for the rjb rates and the last three lines 
for T. The T(5) is a well-studied particle, and so we can use the measured decay rates 
to estimate ag and clb s - While the total width of the T(6) is known to within 20%, the 
individual hadronic channels of interest are currently unknown. Therefore, we assume its 
branching ratios to the various B and B s final states are the same as those of the T(5), 
which should not induce a significant further error in our admittedly approximate treatment 
of the threshold region. 

To calculate the a decay rate into the states of Eq. \A1\ we replace t#, with a, and the use 
the effective coupling on the RHS, 

\a Ba \ 2 = |ob5| 2 |c 5 | 2 + |ob6| 2 |c 6 | 2 (A3) 

I 1 2 I |2||2 I |2||2 

\ a B s a\ = \0-B 3 5\ \C5\ + \dB 3 6\ | Ce | , 

where c§ t 6 are given Eq. [2HJ We have made an assumption about the relative phase between 
aB5 and a Be, but our final answer will have a small dependence on this as typically only 
one of c 5 or c 6 is large. We plot the inclusive decay rate of a to bb in Fig. El seeing two 
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FIG. 6: Decay rates in GeV of the a to bb inclusive. (Blue): Partonic QCD computation assuming 
mb = m B = 5.279 GeV with NLO correction given in Ref. [44|. (Red): Decay rate of a to 6-flavored 
mesons from its mixing with % states. (Yellow): Interpolation given by Eq. lA4l to cover region near 
B meson threshold. (Green): Mixing calculation without B s decays, which violate the assumption 
of uniform coupling, as a in Eq. IA1I at the order of magnitude level. 



pronounced peaks. We use the interpolation function of [44| to take us from the decays due 
to mixing near threshold to the partonic calculation at larger m a : 



T(a — > bb) = T(a — > B, B s mesons via rj b ) + r part . 



onic-NLO 



[a — > bb) ( 1 — exp 



a 



(A4) 

where x = \Jl — — B+ J^ B *" > , b = 5 so that the mixing contribution dominates just above 
threshold, and a = 0.32 is a sum-rule inspired normalization factor so that the integral of 
the interpolating function above threshold equals that of NLO partonic calculation. 
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